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Abstract 

In deep-inelastic processes the heavy flavor Wilson coefficients factorize for Q 2 m? into 
the light flavor Wilson coefficients of the corresponding process and the massive operator 
matrix elements (OMEs). We calculate the 0(a 2 ) and 0(af) massive OME for the flavor 
non-singlet transversity distribution. At 0(a 2 ) the OME is obtained for general values 
of the Mellin variable N, while at O(a^) the moments N = 1 to 13 are computed. The 
terms oc Tp of the 3-loop transversity anomalous dimension are obtained and results in the 
literature are confirmed. We discuss the relation of these contributions to the Soffer bound 
for transversity. 



1 Introduction 



The transversity distribution At/(x, Q 2 ) is one of the three possible quarkonic twist- 2 parton 
distributions besides the unpolarized quark density f(x,Q 2 ) and the longitudinally polarized 
density Af(x,Q 2 ). Unlike the latter distributions, it cannot be measured in inclusive deeply 
inelastic scattering since the corresponding contribution is oc m 2 /Q 2 , [1], with m q a light quark 
mass and Q 2 the virtuality of the exchanged gauge boson. It can be extracted from deep- 
inelastic scattering studying isolated meson production, also called semi- inclusive deeply-inelastic 
scattering (SIDIS), [2,3], and in the polarized Drell-Yan process [3-5]. 1 Measurements of the 
transversity distribution in different polarized hard scattering processes are currently performed 
or in preparation [7] . In the past phenomenological models for the transversity distribution were 
developed based on bag-like models, chiral models, light-cone models, spectator models, and 
non-perturbative QCD calculations, cf. [8] . The main behaviour of the distributions is that they 
vanish by some power law both at small and large values of Bjorken x and exhibit a shifted 
bell-like shape. First attempts to extract the distributions out of data were made in Refs. [9]. 
The moments of the transversity distribution can be measured in lattice simulations, which help 
to constrain it ab initio. First results were given in Refs. [10]. From these investigations there is 
evidence, that the up-quark distribution is positive while the down-quark distribution is negative, 
with first moments between 0.85 ... 1.0 and -0.20 . . . -0.24, respectively. 

The scaling violations of the transversity distribution were explored in leading-, [5,11-13] 2 , 
and next-to-leading order, [15-17]. 3 In Ref. [12] also the method proposed in [19] was used 
to calculate the anomalous dimension. At three-loop order the moments N = 1 to 8 for the 
anomalous dimension are known [20]. For the calculation of the scattering cross sections also 
the corresponding Wilson coefficients have to be known. In case of SIDIS these corrections have 
not yet been calculated. For the transversely polarized Drell-Yan process the 0(a s ) Wilson 
coefficient was derived in Ref. [17] based on [21] and at higher orders the contributions due to 
soft-resummation are available [22]. 

The scattering cross sections dominated by the transversity distribution receive heavy fla- 
vor corrections, although transversity itself is a flavor non-singlet distribution. These contri- 
butions reside in the corresponding Wilson coefficients. In deep-inelastic processes the heavy 
flavor Wilson coefficients factorize into massive operator matrix elements (OMEs) and the 
light flavor Wilson coefficients at large enough momentum transfer Q 2 3> m 2 , as was shown 
in Ref. [23], with m the heavy quark mass. In this way all contributions except the power 
corrections (m 2 /Q 2 ) k , k > 1 can be calculated. The massive OMEs derive from the twist 2 
operators emerging in the light-cone expansion between on-shell states and are process inde- 
pendent quantities. The formalism proposed in Ref. [23] has been applied successfully to calculate 
the asymptotic heavy flavor Wilson coefficients at 0(a 2 s ) [23-25] in unpolarized and polarized 
deep-inelastic scattering. For Fl(x,Q 2 ) the asymptotic heavy flavor corrections to 0(a 3 ) were 
calculated in Ref. [26]. A series of Mellin moments for the asymptotic heavy flavor Wilson coef- 
ficients contributing to the structure function F 2 (x,Q 2 ) at 0(a^) have recently been computed 
in Refs. [27,28]. 

In the present paper we apply this formalism to the tensor-operator defining the flavor non- 
singlet transversity distribution, and limit the consideration to contributions of twist 2 and the 
collinear parton model. We calculate the 0(a 2 ) corrections for the flavor non-singlet OME of 

1 For a review see Ref. [6]. 

2 The small x limit of the LO anomalous dimension was calculated in [14]. 
3 For calculations in the non-forward case see [12, 18]. 
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transversity for general values of the Mellin variable N. At 0{a*) the OME is computed for 
individual Mellin moments N — 1 to 13. The 2-loop calculation verifies the TV-terms of the 
transversity anomalous dimension of former NLO calculations [15-17]. In the 3-loop calculation 
we obtain the moments for the complete 2-loop anomalous dimension, which appears in the 
double pole term in the dimensional parameter e = D — 4. Furthermore, the TV-contributions to 
the 3-loop anomalous dimension are obtained from the single pole term, which can be compared 
to the results in [20] for N = 1 to 8, while the TV-terms of the anomalous dimension for N = 9 
to 13 are new. The results for the massive OME for transversity given in the present paper are 
related to future lattice simulations with (2+1+1)-, resp. (2+l)-dynamical fermions. The heavy 
flavor contributions are also of relevance for the Soffer bound for transversity [29]. 

The paper is organized as follows. In Section 2 we summarize the main relations for semi- 
inclusive scattering cross sections in the leading twist approximation from which the transversity 
distribution can be determined. Here, as in the case of inclusive deep-inelastic scattering, tagging 
on charm-mesons allows to measure the charm contribution directly in high-luminosity experi- 
ments. The method to calculate the heavy flavor corrections in the asymptotic region is briefly 
described. In Section 3 we calculate the O(otf) massive operator matrix element. The Mellin 
moments of the OME at 0(a 3 s ) are computed in Section 4. In Section 5 we discuss the heavy 
flavor contributions to the Soffer bound and Section 6 contains the conclusions. In the Appendix 
we summarize the Tp-parts of the 3-loop anomalous dimension for transversity and the moments 
of the constant part 0(e) of the un-renormalized 0(<X?) massive OME for the Mellin moments 
N — 1 to 13. For details concerning the calculation and renormalization of massive non-singlet 
OMEs we refer to Ref. [27]. 



2 Basic Formalism 

The transversity distribution 

A T f(x,Q 2 ) = f(x,Q 2 )-f^(x,Q 2 ) (1) 

contributes to a large variety of scattering processes, cf. [6]. Here j (j) denote the transverse spin 
directions. Eq. (1) describes the transversity distribution obtained in the light-cone expansion 
at twist 2 or in the collinear parton model. For other phenomenological applications one may 
introduce /cj_-effects for this distribution, [6]. This, however, has consequences for the twist 
expansion and the renormalization of the corresponding processes, when calculating them to 
higher orders. We will therefore restrict the analysis to the level of twist 2 and consider only 
processes which are free of /cj_-effects, or after these were integrated out in the final state. 

For semi-inclusive deeply inelastic charged lepton-nucleon scattering IN — > l'h + X the Born 
cross section, after the P ^-integration, is given by, [6], 

d3(T £ [1 + (1 -Vf] F a (x,Q 2 )D a (z,Q 2 ) 



dxdydz Q 4 a 2 



a=q,q 



-(1 - sOlSjSfcJ cos (0 5 + 5 J A T F a (x, Q 2 )A T D a {z, Q 2 ) j . (2) 

Here, in addition to the Bjorken variables x and y, the fragmentation variable z occurs. Sj_ and 
S/ij_ are the transverse spin vectors of the incoming nucleon N and the measured hadron h. The 
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angles <ps,s h are measured in the plane perpendicular to the 7*iV (z-) axis between the x-axis 
and the respective vector. The transversity distribution can be obtained from Eq. (2) for a 
transversely polarized hadron h by measuring its polarization. The functions Fi,D h A T Fj, A T Di 
are given by 

F t (x,Q 2 ) = CifoQ^O/ifoQ 2 ) , (3) 
D, t (z,Q 2 ) = C t (z,Q 2 )® D^z.Q 2 ) , (4) 
A T F,(x,g 2 ) = A T C t (x,Q 2 )® A T fi(x,Q 2 ) , (5) 
A T L>i(z, Q 2 ) = A T C t (z,Q 2 )®A T D t (z,Q 2 ). (6) 

Here, <8> denotes the Mellin convolution, Di, A^Di are the fragmentation functions and 
Cj, Cj, A T Cj, A T Ci are the corresponding space- and time-like Wilson coefficients. The Wil- 
son coefficient for transversity, A T Ci(x, Q 2 ), contains light (A T Cj) and heavy flavor (A T Hj) 
contributions 

AtC^x, Q 2 ) = A T Ci(x, Q 2 ) + A T H t (x, Q 2 ) . (7) 

For brevity we dropped arguments like m 2 , the factorization scale, /z 2 , and the number of light 
flavors, Nf, in Eq. (7). 

Eq. (2) holds for spin-1/2 hadrons in the final state, but the transversity distribution may 
also be measured in the lepto-production process of spin-1 hadrons, [30]. In this case, the 
Pfcj_-integrated Born cross section reads 

( '° l?:a >inU,. • o, (r ) S. .S, r (1 //) ^ (7 VA, /•;(,-. g 2 )//,,,.,^:.^ 2 ) . (8) 



dxdydz xyQ 2 



i=q,q 



Here, the polarization state of a spin-1 particle is described by a tensor with five independent 
components, [31]. 4>lt denotes the azimuthal angle of Slt, with 

\Slt\ = yJ(S* LT ) 2 + (Sl T f . (9) 

H a ,i,LT(z,Q 2 ) is a T- and chirally odd twist-2 fragmentation function at vanishing k±. Process 
(8) has the advantage that the transverse polarization of the produced hadron can be measured 
from its decay products. 

The transversity distribution can also be measured in the transversely polarized Drell-Yan 
process using the polarization asymmetry, see Refs. [17,21,22]. However, the SIDIS processes 
have the advantage that in high luminosity experiments, cf. [32], the heavy flavor contributions 
can be tagged like in deep-inelastic scattering. This is not the case for the Drell-Yan process, 
where the heavy flavor effects appear as inclusive radiative corrections in the Wilson coefficients. 
We will therefore mainly consider SIDIS in the following. 

As was shown in Ref. [23], in the region Q 2 ^> m 2 all non-power contributions to the heavy 
quark Wilson coefficients obey factorization relations. In the general flavor non-singlet case one 
obtains for Nf light and one heavy quark 

flr™- (*. f ) = (x, J, if, + 1) 45, (*, J) - cj» J, iv,) , no) 
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where is a light flavor Wilson coefficient and A^q is the corresponding massive operator 
matrix element, cf. [23,25,27], with 



NS 



k=l 



m 



/// 



(g|0 NS |g) = 5(l- a ;) + ^^(^)<! Q NS 



k=2 



m 



(11) 

(12) 



Here a s (fi 2 ) = a s (fi 2 ) / (An) denotes the strong coupling constant and \q) are light quark states, 
with on-shell momenta. The local flavor non-singlet twist- 2 operator for transversity is given by 



qTR,NS,h,hi,... ,hn ( z j 



;N-1 



s 



q(z)a^D^ ...D» N —q(z 



— Trace Terms , 



(13) 



with a*™ = (if 2) [Yl u -YY], K the Gell-Mann matrices for S£/(3) flavor , D v the covariant 
derivative in QCD, q(q) denote the quark and antiquark fields, and the operator S symmetrizes 
the Lorentz indices. Note that in Eq. (10) the heavy quark degrees of freedom are all contained 
in the process independent OMEs. 

In case of transversity one obtains the following representation for the heavy flavor Wilson 
coefficient 4 after expanding Eq. (10) up to 0(a^) 

A T H^(N f + 1) = a 2 (N f + 1) [A T A^ S + A T Cf (A^ 

+ a 3 a (N f + 1) \A T A^ s (N f + 1) + A T ^ Q NS A T Cf 



+A T Cf (iV» 

Here we made the Aj-dependence explicit and use the notation 

/(TV,) = f(N f + 1) - f(N f ) . 



(14) 



(15) 



We dropped all arguments like x,N,m 2 ,fi 2 , which are understood implicitly. Additionally, 
Eq. (14) is written in Mellin space, in which we will work from now on, if not stated other- 
wise. The assignment of the differing arguments in Nf in Eq. (14) is necessary to project onto 
the heavy quark part. 

Following Ref. [27] we consider the Green's function G^J : Q ,m which is obtained by contracting 
the matrix element of the local operator (13) with the source term J N = A M . . . A MJV 



u( P ,s)G^ m X r u(p,s) = J N (qi(p) I O'q^l 



TFLNS 



,111,... ,HN 



q J (p)) 



Q ? 



(16) 



where p and s denote the 4-vectors of the momentum and spin of the external light quark line, 
u(p, s) is the corresponding bi-spinor, A. A = 0, and Q labels the heavy quark contribution. The 
un-renormalized Green's function has the following Lorentz structure 



Aij,TR,NS 



= * y (A • p) N - x \A p o^A T A qqtQ (— , e, Nj + Cl A^ + c 2 ^ + <#fi> 

+c 4 Aj^ + c 5 A i>p^ , (17) 



4 Apparently, the light flavor Wilson coefficients for SIDIS were not yet calculated even at 0(a s ), although 
this calculation and the corresponding soft-exponentiation should be straightforward. 
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with unphysical constants Cfc|fc=i...5 and fh the un-renormalized heavy quark mass. The un- 
renormalized massive OME is then obtained in Mellin space via the projection 

NS A 2 £ij 

^A qq , Q {^N) = ~i mA J N+1{D _ 2) {Tr[A i> P^ m ] - A.pTr^G^ 5 ] 

+.A. P Tr[ w "G^ T J NS ]} . (18) 

Here N c denotes the number of colors. For the renormalization procedure and different steps to 
the final representation of the massive OME in the MS-scheme we refer to Ref. [27]. Note that 
the renormalization of the heavy quark mass is carried out in the on-mass-shell scheme. 



3 The 0(a 2 s ) Massive Operator Matrix Element 

After mass renormalization the massive flavor non-singlet OME for transversity at 0(a 2 s ) is given 
by [23,25] 



S 



in 



/// 



2\ 6 



(1),TR 



(2),TR _(2),TR 
+ a qq,Q + £a qq,Q 



(19) 



Here we dropped all arguments on the r.h.s. S £ is the spherical factor which occurs due to 



dimensional regularization and is set to one in the MS-scheme. 



(fc),TR 



111 



(N) denote the (k + 1)- 



loop anomalous dimensions for the non-singlet composite operator (13). Note that as in Ref. [27] 
we define the anomalous dimension corresponding to an operator Z-factor via 



d 

7 = n-^hi(Z{n)) 



(20) 



/?o,q denotes the heavy flavor contribution to the /9-function in lowest order, 
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(21) 



with Tp = 1/2. Eq. (19) has been expanded up to 0(e) since the coefficient a qq n enters the 
3-loop OME via renormalization. The renormalized OME is given in Mellin space by 



A 4 (2),NS,MS 
A T A qq,Q 



(N) 



n (0),TR 2 

4 KyP) 



-(1),TR 2 



+ a 



(2),TR 
qq,Q 



(0),TR 



<2 (22) 



in the MS-scheme and k > 2, k e N denotes the Riemann (^-function at integer arguments. 
The calculation of the 2-loop OME in terms of Feynman-parameter integrals is straightforward, 
see [25]. For the anomalous dimensions 7g°' ) ' TR and 7^' TR we obtain 



2C F [-3 + 4^] , 



a(i),tr 

Iqq 



(N) = fc F T F 



3S 2 - 5S 1 + - 



(23) 
(24) 



with C F = (N 2 — 1)/(2N C ), confirming earlier results, [15-17]. Here, Sk = Sk{N) denote the 
single harmonic sums, [33]. The finite and 0(e) contributions of the un-renormalized OME, 
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Eq. (19), read 



a 



(2),TR 
qq,Q 



(N) 



C F T, 



pip' 



40 



5, 



224 8 ' 



Si + 2C 2 



(24 + 73 N + 73 N 2 
187V (N+ 1) 



_(2),TR 
a qq,Q 



(N) 



C F T F { - 



656 20 8 
"sT + Y C2+ 9 Cs 



5*1 + 



112 4 
^7" + 3 G 



9 20 9 



4 1 2 (-144-48iV + 757iV 2 + 1034iV 3 + 517iV 4 ) 
+ 3^ 4+ 6 C2 + 3 C3+ 216iV 2 (iV+l) 2 



The renormalized 2-loop massive OME (22) then becomes 



(2),NS,MS 

qq,Q 



(N) = C F T F < 



-3 S ' + 2 



m 



¥ Sl + 3 + T s ' 2 



8 40 224 24 + 73iV + 73iV 2 
3 3 + 9 2 " ^7" 61 + 18iV(iV + l) 



(25) 



(26) 



(27) 



Corresponding quantities for vector currents were calculated in Refs. [23,25]. In the limit N — > oo 
7 S ) (iV),7g ) (iV),a^ Q (Ar),a^ Q (Ar) and A® Q (N) in the vector and transversity case approach 
each other. This has also been observed for the 2-loop transversity anomalous dimension in 
Ref. [15]. 



4 The O(a^) Massive Operator Matrix Element 

The renormalized OME for transversity at 0(a^) has the same structure as the flavor non-singlet 
OME in the case of vector currents, [27]. In Mellin space it is given by 



A 4 (3),NS,MS 
^\q,Q 



(AO = 



(0),TR fl 
Iqq A),Q 



. (m 



P + 2p 0tQ )\n^— )+ 7 ^2 7 5i MK A», 



4 | " '<"/ 

2. 



>,Q 



-2^{Po + P , Q ) + A, Q 7S?' TR | ln 2 (^) + ±|t£>.™ 



+4oSq TR (A + A),g) - 7^ 
+2^/30^?'™ + ^S 0) 7«' TR + 2^S- 1 )4J Q TR + ^ (28) 



o/ (0),TR^ ^ , (l),TRo , 

(2) 7gg PoPo,Q(,3 7gg Po,q(,2 
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in the MS-scheme, performing mass renormalization in the on- mass-shell scheme. Here the 
expansion coefficients of the /^-function and the mass renormalization constants are, cf. [27,34,35], 

(29) 

(30) 

(31) 

(32) 

(33) 
(34) 

(35) 



Po = 


jC A - \t f n, , 


^ - 
Pi,Q — 


— 4 ( -Ha 4- CtA 7V. 
^ 1 <^ A ' ^ F J F 1 


8 {1) - 


32 
y 


B (2) - 


86 31 


5m[ 1 - > = 




Sm^ = 


-4C F , 


5m^ = 





(\) TR 

with Ca = N c , and the NLO anomalous dimension 7^ reads, cf. [15-17], 

C 2 F (AS 2 - 8S 1 - 1) 
C A 



+8C F C F - 



AS±S 2 - 8S 1 S^ 2 + Si- AS 3 - 4S_ 3 + ^S 2 



+8S. 



2,1 



■1) 



N 



N(N + 1) 



+C F C A ( -165x5 2 - fs 2 + ^ - y 

, ,'32 „ 160 „ 4 

+C F T F N f — S 2 5i + - 

'3 9 3 



(36) 



All contributions to Eq. (28) are known for general values of N, except of 7gg ) ' TR and o^q R , 
the constant contribution to the un-renormalized 3-loop massive OME. Similarly to the vector 
case, Ref. [27], we calculate A t A^q s ' MS (jV) for a fixed number of Mellin moments. The Feynman 
diagrams were generated by a code based on QGRAF [36] and the color algebra was performed 
using color, h [37]. The computation is based on FORM [38] codes using MAT AD [39]. Since the 
projector in Eq. (18) has to be applied we can calculate the moments N — 1 to 13, i.e. one 
moment less than in the vector case in Ref. [27], given the complexity of the problem and the 
computer resources presently available. The computation time amounted to about 9 days. The 

(2) TR 

contributions to the 3-loop transversity anomalous dimension, y qq (N), and the constant part 

( 3) TR 

of the un-renormalized massive transversity OME a qq \Q (N) are given in the Appendix. In 

Figure 1 the numerical values of a> q 3 q Q R (N) are compared to those of o^q(N), Ref. [27]. As has 
been observed at O(o%) already, for larger values of N both quantities approach each other at 
0{a 3 s ). This also applies to 7 g ),(V ' TR) (N). 

The present calculation confirms the Tp-parts of the transversity 3-loop anomalous dimension, 
which was calculated for N = 1 to 8 in Refs. [20] for the first time. We also present the 
moments N = 9 to 13. As a by-product of the present calculation the complete NLO anomalous 
dimension [15-17] is confirmed for the moments N — 1 to 13. 
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Finally, we show as examples the first moments of the MS-renormalized 0(a^) massive 
transversity OME. Unlike the case for the vector current, the first moment does not vanish, 
since there is no conservation law to enforce this. One obtains 



(3),NS,MS 



T/i qq,Q (!) 
104 



27 

496 
~81 



^{(lic,-15 W + 2 ))^(f) + (fc-^ 



In 



m 



16 „ 278, 

— B 4 + 24( 4 Cs + 

3 s 9 81 



27 

(3),NS,MS 
<?<?,<2 



Cs- 
(2) = 



/44 
\~9 



C I { i-y/,( i \/-2)) 



-87>) ln 2 (^) + [(l5 + 48C 3 )C F + (-^ - 48C 3 )CU - ™N f T F 




13/ 34135\ / 6556 128. \_„ /2746 224 
729 ) A+ \ "729~ + ~2T'V~ X ~ 1 ' V 729 



(37) 



V 



2 ) + (-> 



196 



+ 



496 

~2T 

533 



In 



(™ ) + (-16H, + nU - + ^§)C F+ (8B 4 - 72 C4 
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C 3 -56)C A + (- 



1988 128 \ 
^T + T~ C V TFiV/ ' V27 



27 

/338 _ 224 



(38) 



with 



-4( 2 ln 2 (2) + | ln 4 (2) - ^4 + 16Li 4 Q 



(39) 



-1.7628000871 



The structure of the massive OME is similar to the result for the unpolarized case, cf. Ref. [27], 
Eq. (5.57). We checked the moments N — 1 to 4 keeping the complete dependence on the gauge- 
parameter £ and find that it cancels in the final result. We observe that the massive OME do 
not depend on (2 5 , as is also the case for the various massive OMEs which were calculated for 
vector currents in Ref. [27]. The results for the massive OME for the moments N = 1 to 13 and 
the quantities listed in the Appendix are attached to this paper in FORM- format. 

Since the light flavor Wilson coefficients for the processes from which the transversity dis- 
tribution can be extracted are not known to 2- and 3-loop order, phenomenological studies 
on the effect of the heavy flavor contributions cannot yet be performed. However, our results 
can be used in comparisons with upcoming lattice simulations of operator matrix elements with 
(2+l+l)-dynamical fermions including the charm quark. 



5 Thc combination of multiple zeta values B4 is characteristic for quantities depending on a single mass scale. 
In this specific combination (Rvalues at even integer argument contribute. 
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5 Remarks on the Soffer Bound 



If the Soffer inequality [29] 



\A T f(x,Ql)\<-[f(x,Ql) + Af(x,Ql)] 



(40) 



holds for the non-perturbative parton distribution functions at a given scale Q% one may check 
its generalization at the level of the corresponding structure functions. In the light-flavor case 
this has been investigated to 0(a s ) for the Drell-Yan process in Ref. [17]. For the heavy flavor 
corrections studied in the present paper one investigates 



\A T F(x,Q 2 )\ < i [F(x,Q 2 ) + AF(x,Q 2 )} , 



(41) 



where the structure functions are given in Eqs. (3,5,7) and by corresponding relations in the 
longitudinally polarized case. One may try to separate the evolution effects in the parton distri- 
bution functions from those of the Wilson coefficients. 

The solution of the non-singlet evolution equation for the parton distribution f NS (N, Q 2 ) in 
Mellin space for Nj massless flavors reads to 3- loop order, cf . [40] , 



f m (N,Q 2 ) 



E{N,Q\Ql) f^(N,Q 2 )=[- 



(0),NS,. n , a 



a 
a 

1 



1 - ^-(a - a ) 
Po 



"7, 



E(N,Q 2 ,Q 2 ) f NS (N,Ql) 
Po 



2,9o 
1 



(a 2 - al) 



7jf NS (iV) + |7«' NS (iV) - (| - |) l^W 



W.NS/yvrX _ fa (0),NS 



[a - a ) 2 




f NS (N,Ql) 



(42) 



where a = a{Ql) and 
34 



A = "" 2 
2857 



20 



3 C A - AC F T F N f - —C A T F N f 



(43) 



ft = 



54 



C\ + 2C F TpNf 



205 



1415 



44 



158 



C' A T F N f + —CpTpNj + —C A T^N: 



(44) 



cf. [41]. The moments of the anomalous dimensions for vector currents are given in Refs. [42]. 
The evolution operator in the unpolarized and the longitudinally polarized case are the same 
due to a Ward identity. Therefore it is sufficient to investigate the relation 



\A T E(N,Q 2 )\ <E V (N,Q 2 ) 



(45) 



Up to the 0(a s ) corrections (NLO) the validity of this inequality was shown in 
[17]. Beyond this level only a finite number of Mellin moments can be compared for 
E TR (N, Q 2 , Ql)/E Y (N : Q 2 , Ql), for which the 3-loop transversity anomalous dimension is known 
[20], expanding up to 0(a 2 ). This quantity is shown in Figure 2 for the 2- and 3-loop case. The 
corresponding correction preserves the Soffer bound for characteristic values of a s . 
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Turning to the effect of the heavy flavor Wilson coefficient in the asymptotic region, Eq. (14), 
we have to limit the investigation to the massive operator matrix elements since the corresponding 
light flavor Wilson coefficients were not yet calculated. In Figure 3 we show the difference 



A (2),NS,MS / \ A . W ,lNO,M3/ N 

A qq,Q K X ) ^T* qqtQ [X) 



.(2),NS,MS, 



C F T F {l-x){ ^ln 2 



m 



/// 



In 



m 



/// 



+ - 



22 



ln 2 (a;) + — \n(x) + 



116 
IT 



(46) 



for a series of values Q 2 = /i 2 . At large scales A q 2 q Q S ' MS (x) — A T A q 2 q Q S ' MS (x) is positive and 
descending towards x ~ 1, while at lower scales also negative values are reached in the interme- 
diate region of x. The difference is always positive in the small x region. To maintain the Soffer 
bound the light flavor Wilson coefficients have to compensate the negative contributions. In 



,(2),NS,MS A a (2),NS,MS 
\q,Q _-^TA qqQ 

The behaviour of A q f^' US - A T A^ MS , Fi gure 5, is quite similar to that shown in Figure 4 



Figure 4 we show A qq Q ' lv '° — A T A { qq Q D '' v '° in Mellin space, where also a sign change is obtained. 



and a corresponding behaviour of A 



(3),NS,MS l 
qq,Q 



X 



A T A 



(3),NS,MS, 



x) to the one found at 0(a 2 ] 



is 



expected. From the knowledge of the massive OMEs alone a conclusion on the validity of the 
Soffer bound at the level of the structure functions can not be drawn before the light flavor 
Wilson coefficients have been computed. 



6 Conclusions 

We calculated the flavor non-singlet massive OME for transversity at 2-loop order and for the 
Mellin moments TV = 1 to 13 at 3-loop order. For large scales Q 2 ^> m 2 the heavy flavor 
Wilson coefficient can be determined from the light flavor Wilson coefficients and the respective 
process independent massive operator matrix element computed in the present paper. For flavor 
non-singlet quantities the heavy flavor corrections start at 0(a 2 ). The measurement of the 
corresponding scattering cross sections requires high luminosity. In the present calculation we 
have verified the Tp-parts of the 3-loop transversity anomalous dimension for the moments N — 1 
to 8 and extended this part up to N = 13. As a general observation we found that both the 
anomalous dimension and the expansion coefficients in e computed in the present calculation 
for transversity approach those in the vector case for large values of the Mellin parameter N. 
We investigated the compatibility of the results of the present calculation with the Soffer bound 
on the level of structure functions. While for the evolution operator the Soffer bound is obeyed 
to 3-loop order, a final conclusion cannot be drawn for the massive operator matrix element at 
0(a 2 ) and 0(a^) alone concerning the massive Wilson coefficients for the whole phase space, 
due to a sign change for ^^q 15 — A T A^ q 'Q S at lower scales of Q 2 and medium values of x. A 
firm conclusion can only be drawn after the yet unknown massless Wilson coefficients have been 
computed. 
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7 Appendix 

The Tp-contributions to the 3-loop anomalous dimensions for N = 1 to 13 
7^ TR (l) = Q 



^fTf 



are given by : 

,3 



7S } ' TR (2) = C F T F 
7^ TR (3) = C F T F 
7S } ' TR (4) = C F T F 



7jf TR (5) = C F T F 
2944 



-±T F (2N f + \) -^c A +™C F + 32(C F -C A )( l 

T F (2N f + l)-^C A -mC F + 96(C F - C A )( 3 

'-™*T F (2N, + 1) -^C A - 2 -^C F + ^(C F - C A )( 3 
243 V ' ; 243 243 3 V A 

14 722 , , 199723 66443 512 

"1215" Tp{2Nf + ^430~ ° A " "486" ° F + ~ {Cp ~ ° a)0 

-^T F (2N f + l)- 5 -^C A -^^C F 
30375 FK i ' 60750 A 303750 F 



+ 



7S ),TR (6) = G 



7 F T F 
1088 



3209758 m/ , T , 3682664 ^ 18622301 

T F (2N f + 1) CU i 

212625 V ' ; 42525 101250 



-r - (C F — C A )( 3 





7S ),TR (7) = C fTi 



F 1 F 
8256 



168501142^ /oAr 1844723441 49282560541 

i f ( 2./V f + 1 ) G4 C f 

10418625 v 1 ' 20837250 243101250 



+ ^(Cf-O00 

7^' TR (8) = CfT, 
8816 



711801943 6056338297 849420853541 

41674500 F ^ /+ ' 66679200 A ~ 3889620000 



(c F - c A )(, 3 



(47) 
(48) 
(49) 

(50) 

(51) 



(52) 



(53) 



7£>' ra (9) = CM 



F 



7^' TR (10) = C fTi 



7 



F 

87856 



(54) 

20096458061 119131812533 24479706761047 

" 1125211500 F ^ f+ ' ~ 1285956000 A ~ 105019740000 F 

(55) 

229508848783 4264058299021 25800817445759 

" 12377326500 ; + ^ ~ 45008460000 A ~ 105019740000 F 



(56) 
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if TR (n) = c fTi 



F±F 



28677274464343 /nAT 75010870835743^ 

■ l F (2J\f + 1) ^or^ooor^r^ ^ A 



1497656506500 



778003380000 



a(2),tr 

Iqq 



(12) 



A(2),TR 
Iqq 



396383896707569599 ^ 1006736,,, 

Of H ( O f — G 4 ) Cq 

1537594013340000 3465 v ; 

383379490933459 38283693844132279 
19469534584500 F ^ /+ ^~ 389390691690000 A 

1237841854306528417 1043696 
4612782040020000 F + 3465 ^ F ~ A ^ 3 



CpTp 



(13) = CpTp 



3465 

66409807459266571 6571493644375020121 
3290351344780500 F ^ /+ ' 65807026895610000 A 

36713319015407141570017 14011568 
131745667845011220000 F + 45045 ^ F ~ A ^ 3 



13 



The constant parts a gq Q (N) of the massive 3-loop OME for iV = 1 to 13 are given 



a 



(3),TR 
qq,Q 



1) 



a 



(3),TR 
qq,Q 



(2) 



(3),TR 

<?<?><? 



a 



(3),TR 
qq,Q 



(4) 



a 



(3),TR 
qq,Q 



(5) 



CpT F 



Asl t+ S B 2At 61 r 26441 V 



52 112 
"27 + ~2T 



+ 



+ 



104 

~2T 

278 



6548 
729~ 



15850\ 
729 J 



256 . 



/ -IF 



, 49 . 15715 16 „ c\ a /- 1 - , 

— Cs + — C2 + #4 + 24 C 4 Cj? 

9 3 162 3 1 



CpTp 



577 , n „ ^ 1 1043 

Cs + 8 5 4 - 72 G + - C2 + 

9 s s 3 162 



112 _ 

[ -4C 2 + — Cs 



4390^ 



+ 



+ 



-8C2 
310 



81 

1388 _ 256 
~81 9~^ 3 



C3 + 33C2 + 



10255 



(3) = C F T F 



3 

/ 40001 



54 



N f T F 

16 B 4 + 72 C 4 ) C F 



405 



, 104 r, 121. 327967\ ^ 

C 3 + — 5 4 - 104 C 4 + — C2 + ^tt^t ) CU 



9 



81 



21870 J 



+ 
+ 

+ 
C F T F 



452 1456 168704 \ 



81 
904 



81 

1354 
"9" 



2187 



. 3821 
Cs + -5^" C2 + 



81 

1170943 208 



81 



4374 



9 



5 4 + 104C 4 C F 



752112 



405 



Cs + 



128 



250 



4400353 \ 



9 ^-128C 4 + — C2+ 21g700 



+ 
+ 

+ 



554 _ 1792 _ 20731907\ Ar m 
C2 + — Cs - 010 _ J NfTp 



81 
1108 



81 
556 



C2 



81 *° 218700 J 
3195707 4096 



109350 



81 



Cs T F 



Cs + ^C 2 + ^^-^^ + 128C 4 )C F 



81 



7442628 



736 



174960 



736 



9 



8488 



1436867309 \ 



2835 + 45 Bi 5 ^ + 2025 ^ + 765 15000 



■ 



+ 



15962 A 10304 A 596707139\ 

■ C2 + Cs - 



2025 



405 



5467500 



) 



NfTp 
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a 



(3),TR 
qq,Q 



(6) 



a 



(3),TR 
qq,Q 



(7) 



(3),TR 
l qq,Q 



a 



(3),TR 
qq,Q 



(9) 



31924 92220539 23552 



+ 



2025 ' 2733750 

47932 A 662674 A 

C3 H C2 + 

225 S 10125 S 



405 

40410914719 1472 



109350000 



(64) 



CfTf 



f 172138 . 272 „ 816 . 10837 . 807041747\ ^ 

C 3 + — B A - — u + C2 + roroi rnn j c A 



+ 
+ 

+ 

+ 
+ 



V 945 
17762 



15 



2025 



53581500 



2025 
35524 



2025 
159296 



C2 + 
C2- 



3808 , 32472719011\ m 

G N f T F 

135 Sd 267907500 J } 

5036315611 8704 



133953750 



135 



675 
27982 



. 81181 , 14845987993 544 „ 816 . , 
C 3 + T77T^C2 + „ Arnnnn - ~r^B 4 + — (a )C f 



1125 



36450000 



15 



(65) 



135 
947138 



, 688 „ 6192 A 620686 . 413587780793 \ 
C 3 + — B 4 - — — U + -7777777777- C2 + 



35 



+ 



99225 
1894276 
99225 

8454104 
33075 



u 1376 . 

C2 + —7^ CS 



35 ' 99225 
1727972700289 



52509870000 / 



45 



13127467500 



N f T, 



268946573689 22016 , , 
C2- nr „n„ rn t^t~ C3 T F 



6563733750 



315 



90495089 . 12873570421651 1376 „ 6192 _ , 
Cs + TTV^-r C2 + nn , ne „ rnnnn 7777- B A + — — C 4 C F 



1157625 



(8) = CpTp 



87613 , 2204 „ 
( 3 + —B A 



+ 
+ 

+ 1 



378 

2030251 
" 198450 
2030251 
99225 

9020054 
33075 



105 



29172150000 



6612 > 11372923 , 
C4 + TTT^TT C2 



35 



1587600 



35 35 



91321974347 \ 
112021056000 ) 



(66) 



4408 . 29573247248999 \ 

C2 + -TT^r Cs - n , nMnAannnn ) N f T F 



135 



210039480000 



4618094363399 70528 
C2 — 7777^7^77777777 ttttt - C3 m f 



105019740000 



945 



, 171321401 . 1316283829306051 

C 3 + onronnn C2 + 



2058000 



2800526400000 



4408 6612 
"T05" 54 + ^5" C |( ; 

'9574759 



(67) 



C F T F 



_ 20956 „ 
Cs + B4 



37422 ' 945 
17524721583739067\ 



20956 „ 16154189 „ 
C 4 + nn „ nnn C2 



105 



2041200 



1497161413440000 J 



19369859 , 41912 > 
- I - t^tttttttt; C2 + — r-^Cs 



1786050 



1215 



2534665670688119 
17013197880000 



N f T, 



15 



/ 19369859 397003835114519 670592 \ 
V 893025 C2 ~ 8506598940000 8505 F 

,_ 85698286 131876277049 1013649109952401819 
" ' 297675 C3 + 1500282000 ^ + 2041583745600000 

41912 ^ 20956 



945 105 
CpTp 



261607183 21964 21964 618627019 

C3 + -D4 7773- C4 + =z , , nnnn C2 



935550 945 105 " 71442000 

176834434840947469\ 



7485807067200000 



C 



A 



4072951 , 43928 , 321908083399769663\ _ m 

G H & N f T F 

357210 s 1215 s 2058596943480000 J 1 

4072951 A 50558522757917663 702848 . \ ^ 



178605 " 1029298471740000 8505 

, 3590290 137983320397 11669499797141374121 
' 11907 + 1500282000 + 22457421201600000 

43928 21964 



945 "" 105 

CpTp 



I 3687221539 _ 251684 251684 , 149112401 _ 

C3 + — 777777" B 4 - ——— C 4 + „ nnnnnn C2 



\ 12162150 10395 1155 16038000 

436508000489627050837\ 



11775174516705600000 



1 



/ 514841791 503368 A 40628987857774916423\ Ar _ 

+ & H 6? iVfTp 

V 43222410 S 13365 S 249090230161080000 / ' 

/ 514841791 6396997235105384423 _ 8053888 \ 

" V 21611205 124545115080540000 ~ 93555 J F 

. 452259130, 191230589104127 . 177979311179110818909401 
+ TYTT^T Cs + 777777777^7777777777 C2 + 



1110717 19968753 J2000 ,_ 328799103812625600000 



503368 251684 

-D4 H — , -, rr C4 



10395 1155 
CfT} 



5827712409 , 245210883820358086333 260924 „ 260924 

C2 " A ^ nnnAn ^ A ^ nrnnnn + 77777777- B i ~ „rr < 



S644482000 4783664647411650000 10395 1155 

3971470819 
+ 12162150 ^ 3 ' ' 

/ 7126865031281296825487 521848 , 535118971 „ . , 

+ 1 6 iVf T F 

V 42096248897222520000 13365 s 43222410 s ' 1 



16 



a 



(3),TR 
qq,Q 



( 8349568 535118971 1124652164258976877487 \ 
+ V 93555 ^ 3 ~~ 21611205 ^ 2 ~~ 21048124448611260000 J 

( 260924 2396383721714622551610173 468587596 
+ ^ 1155 ^ 4 + 4274388349564132800000 1440747 ^ 3 

521848 „ 198011292882437 ^ 

B 4 H C2 Cj 

10395 1996875342000 S I 



TV 



(71) 



(13) = C F T F 



/ 15314434459241 
I 1460917458000 



430633219615523278883051 3502892 



6467514603300550800000 



135135 



B, 



3502892 _ 327241423 , 
C 4 + Mrrrn C3 



lev 



+ 



15015 " 935550 

1245167831299024242467303 
7114266063630605880000 



7005784 



93611152819 



+ - 



112092544 



+ 



1216215 
/ 3502892 



93611152819 
3652293645 



173745 " 7304587290 
196897887865971730295303\ 
3557133031815302940000 J 



(2 N f T F 



15015 



C4 + 



70680445585608577308861582893 81735983092 



122080805651901196900800000 



243486243 



7005784 449066258795623169 , ,.. 
-04 H 777777777777 ^o^.^^ C2 I 



135135 



4387135126374000 



(72) 
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Figure 1: The constant part oL q of the un-renormalized flavor non-singlet massive 3-loop OME 
in the vector case [27] and for transversity for Nf = 3. 
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Figure 2: Ratio of the evolution operators .E TR,v (iV), Eq. (42), expanded up to the 0(a s ) terms 
(2 loops) and the 0(a%) terms (3 loops), respectively, as a function of the Mellin variable N, 
with a s o = 0.3. 
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Figure 3: Difference of the massive OMEs A^q S ' MS (x) - A t A ( ^'q S ' MS (x), Eq. (46), for different 
ratios Q 2 /m 2 . 
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Figure 4: Difference of the massive OMEs A®>q S > MS (N) - A t A { ^'q S ' MS (N), Eq. (27), and 
Ref. [23], Eq. (3.35) for different ratios Q 2 /m 2 . 
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Figure 5: Difference of the massive OMEs A { ^'q S ' MS (N) - A t A { ^'q S ' MS (N), Eq. (28), and 
Ref. [27], Eq. (4.17) in the MS-scheme for different ratios Q 2 /m 2 . 
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